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A STATISTICAL TECHNIQUE FOR COMPUTER IDENTIFICATION OF 


OUTLIERS IN MULTIVARIATE DATA 
Ram Swaroop 

Computing and Software, Inc. 
Field Team at Flight Research Center 

and 

William R. Winter 
Flight Research Center 


INTRODUCTION 


The NASA Flight Research Center is engaged in an extensive biomedical research 
and development program. Objectives of this program include advancing the state of 
the art in the medical monitoring of humans in flight (ref. 1); predicting and extending 
the limit of man' s operational capacity in the flight environment; and developing im- 
proved protection, restraint, and life support systems. As a result of this program, 
large quantities of biomedical information are collected in flight, necessitating de- 
pendence on the Flight Research Center's capacity for collecting, reducing, and ana- 
lyzing these data by automatic means. 

Experience has shown that no matter how sophisticated the monitoring, collection, 
and reduction systems, some editing of the biodata is required before they can be ana- 
lyzed statistically. The reduced biodata may contain observations that deviate markedly 
from the rest of the sample. Such observations may be due to errors other than the 
usual random fluctuations characterizing the population to which the data belong, or 
may merely occur too infrequently to be considered in a particular analysis. If, upon 
examination, an observation falls outside a standardized region, it is usually identified 
as an outlier. Outliers often provide useful information. Their identification is im- 
portant not only for improving the analysis but also for indicating anomalies which may 
require further investigation. 

A statistical technique, and the related computer program, for identifying the 
outliers in univariate data was presented in reference 2. A method for identifying 
outliers in multivariate data is derived and demonstrated in this report. This method 
was chosen because of its simplicity and applicability in editing biodata. A program 
for automatic editing was written in FORTRAN IV. Inputs to this program are the 
number of variables, the data set, and the selected level of significance. An example 
is presented to illustrate the use of the method, and a scatter plot of the data is shown. 
The program source listing, user instructions, and a sample output are also presented. 

The program computes and prints the means and standard deviations of all the 
variables before and after the outliers are identified and deleted. A list of the data 
with outliers identified by asterisks is also printed. 


The authors would like to acknowledge the assistance of M. C. Nesel in writing 
the computer program. 


SYMBOLS 


A nonsingular matrix 

^ 05 p,n-p-l cv-level value of F-distribution with p and (n - p - 1) degrees 

’ ’ of freedom 


G 


H/R 

I 


k 

Z 


i 


Np(M,r) 


n 

P 

S 


S.D. 


T 


2 


u,v 

Xj,X. 


X 

Zi 

a 

A. 

1 

A:f; 


normal component of acceleration as experienced by the subject, 

g 

heart rate, beats per minute 
identity matrix 

summation starting from i through k, where i and k are 
integers between 1 and n, and i is less than k 

p-variate normal distribution with mean, p, and covariance 
matrix, r 

sample size 

number of variables 

(p X p) matrix of sums and cross-products of deviations of 
observations from X divided by (n - 1) 

standard deviation 

Hotelling' s statistic 

column vectors of p dimensions 

ith or jth observation vector of p dimensions, where i or j 
ranges from 1 to n 

mean vector computed from n observation vectors 

ith vector obtained by orthogonal transformation of vector Xj 

level of significance 

positive real number corresponding to observation vector Xj, 
where i ranges from 1 to n 

positive real number computed from Fq^p, n - p - 1 data. 

set, to compare with Aj 


2 



random variable related to 


Superscript: 

T transpose 


BRIEF DESCRIPTION OF TECHNIQUE 


Outliers are identified by computing, at the given level of significance, the critical 
value, A^, for the data set and for each observation vector, Xj. If Aj is larger 

than A;j<, observation X- is identified as an outlier. The quantity A^ is a function of 

total sample size, n, number of variables, p, and the F-value for the given level of 
significance, whereas each Aj is a function of the observation Xj and the estimated 

mean and covariance matrix from all the observations. It is assumed that all the ob- 
servations constitute a random sample from a p-variate normal distribution. 


DERIVATION OF TECHNIQUE 


Let X]^, X2, . . . , Xj^ be a random sample of size n from a p-dimensional 

normal distribution, Np(/x, r). The observations will be considered as n greater 

than p + 2 column vectors in a p-dimensional vector space. Consider any (n X n) 
orthogonal matrix, with first two rows as shown, 

1 

1 

yfrUpTHT) 


1 






- 1 ) 


representing the rotation of n-dimensional space so that the observations Xj^, X£, 
..., Xj^ are transformed into vectors Z^, Z2, ..., Zj^, where 



Z 


2 



n 

z 


2y]n{n - 1) 



(Xi - X) 


3 


It may be noted tiiat Zx is distributed as Np(\^/i,r), Z2, Z3, Zjj are all 

distributed as Np(0,r), and all are stochastically independent of one another (ref. 3, 

pp. 50-52). Let S denote the estimate of the covariance matrix F. Then the following 
relation holds: 

n _ _ „ n T 

(n - 1)S = 2 (Xi - X)(Xi - X)T = 2 ZjZj 
1 2 


Define a (p X p) matrix, Sq, such that 

n .p 

(n - 2)So = 2 Z^Z. 

3 


Independence of Z2 and set (Zg, ...» Zj^ implies that Zg is independent of Sq and 
the 


t 2 = Z2\~^Z2 ( 1 ) 


statistic is distributed as Hotelling’ s T^. From the relationship between and F 
(ref. 3, pp. 106-107), it follows that 




- 1 , 


is distributed as 


P(n-2) 

(n - p - 1) P»n - p - 1 

Because Z£ is not independent of S, the preceding distribution does not hold for 

t2 = Z2^S"^Z2 (2) 


and the distribution of must be derived. 
By the preceding definitions 


n 


(n - 1)S — Z2Z2 + 2 ZjZj — Z2Z2 + (n - 2 )Sq 
3 
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or 


(n - 2 )Sq = (n - 1)S - (3) 

To express the relation between and the following lemma is used: 

Lemma; Let A be a (p X p) nonsingular matrix and u, v be p-dimensional vectors. 
Then 


/A a - 1 (A ^uHv'^A 

(A - uv -^ ) = A + ^ ' 


(4) 


1 - v'^A ^u 


Proof; The proof of the lemma is presented in appendix A. 
Applying the result (eq. (4)) of the lemma to equation (3), 


•1 _ n - 2 
~ n - 1 


S-1 + 


S"^Z2Z2’^S“^ 


(n - 1) - Z 2 '^S- 1 z 2 


-1 


Substituting this expression for Sq in equation (1) and applying equation (2), 





n - 2 




S“^ + 


S“^Z2Z2^S ^ 


(n - 1) - Z 2 ^S"^Z 2 




n - 2 
n - 1 


Z2'^S“lZ2 + ^ TA A (Z2'^S"1Z2Z2V1Z2) 

(n - 1 ) - Z2 S ^Z2 


y (5) 


n - 2 
n - 1 


2 

T 


+ 


4 

T 


(n - 1) 


t2 


= - 2)7^ 

(n - 1) - 7^ ^ 

O 

This relation provides the distribution of t , and appropriate probability statements 
can be made. 
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Define 


Ai = (Xj - X)'^S"l(Xj - X) for i = 1, 2, . . . , n 
With no loss in generality, is used. From equation (2) 


and from equation (5) 


2 _ 7 T„-l _ n . 
r - Z 2 S Z 2 - n - 1^1 


= 


(n - 2) 


n 


n - 


(n - 1) 


n 


n - 11 
n(n - 2)Aj^ 

(n - 1)^ - nA-i 


is distributed as 


P(n - 2) 


(n - p - 1) p,n-p-l 

9 

From the distribution of T^, the statement 

n(n - 2)Ai ' 


Probability 


ity - 

I 


(n - 1)2 - nAi “ (n"- P - 1) ^cr,P.n-p-l 


a 


provides criteria for identifying the Z£ (or Xj^) as an outlier at the assigned level of 
significance, a. This statement is equivalent to 


Probability 


^1^ 


p(n 1) 
n(n-P- l)+npF^.p^^_p_l J 


= a 


For significance level a, denote 


^ 1) Fq,.p^j^_p_]^ 

*“n(n-p-l)+npF^.p^^_p_^ 

then Xj^ will be identified as an outlier at a level if 


> A* 
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The quantity (or Zg) was chosen for convenience of the preceding derivation and 
the derivation holds for all Xj. Thus Xj^ will be identified as an outlier at a level of 
significance if A. > A^. 

PROGRAM APPLICATION 


Given a sample of data vectors Xj^, X 2 , • . . , Xjj, the mean vector 

1 

and the estimate of the covariance matrix 

n 

s (Xj - X)(Xi - 

1 

is computed. Then, for assigned a, the critical value 

* n(n - p - 1) + npF^p_„_p_l 

for the data set is computed. Corresponding to each observation vector, Xj, 

Aj = (Xj - X)'^S“^(Xi - X) 

is computed. If Aj > A^, observation vector X. is identified as an outlier at level a. 

The program (appendix B) follows this technique. The output (appendix C) of the 
program contains the data set, Aj, A*, outliers marked by asterisks (*), the number 

of outliers identified, and the level of significance. The output also shows the means 
and standard deviations of the variables before and after the deletion of outliers. The 
required input parameters are: (1) format of the data to be read, (2) number of vari- 
ables, (3) significance level, a, and (4) the data set, formatted as specified. Program 
options allow the user to select either a 5 percent or a 1 percent level of significance 
and to print the names of the variables, if desired. This program is designed so that 
it can be used as a subroutine in other than biodata applications, in engineering and the 
physical sciences, for example. 

The program is particularly useful when large quantities of data are collected and 
the editing must be performed by automatic means. 
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EXAMPLE 


Heart rate, H/R, and normal component of acceleration, G, data from a 66- 
minute flight by a student pilot at the Aerospace Research Pilot School, Edwards Air 
Force Base, Calif. , are used to demonstrate the described technique of computer 
editing of biodata. These data were chosen because centrifuge studies (ref. 4) have 
shown that H/R and G are linearly related. The program was used to identify the 
outliers at a 1 percent level of significance considering H/R and G separately as 
univariate data and together as bivariate data. The computer output for these cases 
is shown in appendix C. 

The results of the two univariate analyses and the one bivariate analysis of the 
same data are presented in figure 1. The point labeled H is identified as an outlier 
on the basis of H/R analysis alone; the point labeled G is identified as an outlier on 
the basis of G alone; and points labeled B are identified as outliers on the basis of 
bivariate analysis of H/R and G. 



cP 


• Outliers 

H Heart rate analysis alone 
G Acceleration analysis alone 
B Heart rate and acceleration bivariate analysis 


1 

2.5 


1 

3.0 


1 _ 

3.5 


.J 

4.0 


Figure 1. Minute heart rate and acceleration data for a 66-minute flight of a student pilot from the Aerospace 
Research Pilot School showing outliers identified by the automatic multivariate outlier technique. 

Bivariate analysis is based on the fact that high H/R is associated with high G, 
whereas univariate analysis cannot take this information into account. For this reason 
the point labeled H was not identified by the bivariate analysis, but was identified as 
an outlier on the basis of imivariate H/R analysis. Also, both points labeled B ap- 
pear not to follow a statistical linear relationship and are identified by the bivariate 
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analysis; however, only one of these points was identified by one of the univariate 
analyses (G alone). This example thus focuses on the fact that the multivariate tech- 
nique, which utilizes the statistical linear relationships between the variables, is 
preferable in identifying outliers in multivariate data. 


CONCLUDING REMARKS 


A statistical technique to identify outliers, or observations which deviate markedly 
from the rest of the sample, in multivariate data at a given level of significance was 
derived. The use of the technique was illustrated by a biodata example. The example 
also demonstrated that the results obtained when each variable was considered sepa- 
rately could be different from the results obtained when the variables were considered 
jointly. The latter technique takes into account the statistical linear relationship be- 
tween the variables and is the preferred method. 

Although this method of detecting and identifying outliers is being used for biodata 
editing at the NASA Flight Research Center, it is also applicable to multivariate data 
encountered in other disciplines, such as engineering and the physical sciences. This 
technique is particularly useful when large quantities of data are collected and the 
editing must be performed by automatic means. 

The program can be used as a subroutine in multivariate analyses. 


I'light Research Center, 

National Aeronautics and Space Administration, 
Edwards, Calif., May 5, 1971. 
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APPENDIX A 


PROOF OF THE LEMMA 


Lemma; If A is a (p X p) nonsingular matrix, and u, v are p-dimensional vectors, 
then 


/A Ta-1 a -1 (A ^u)(v^A 

(A-uv^) ^ = A ^ ^ ^ 

1 - v^A'^u 


Proof; The result is obtained by showing that 

1 /ttT a 


(A - uvT) [a- 1 + (A:iy)(v^A: ) -[ . J 

L 1 - vTa-1u J 


Simplification of the left-hand expression gives 


. . -1 AA ^u(v^A Ta-1 uv^A ^uv^A ^ 

AA ^ + ± : — ' - uv A ^ 


1 - v'^A ^u 


1 - v^A ^u 


or 

I + fuv^A”^ - uv^A ^ + (v"^A~^u)(uv^A“l) - uv^A“^uv^A~^ 

1 - v^A“^u ^ 

T -1 

Because v A u is a scalar, the expression becomes 

1+ ^ — - — [(v^A“^u)(uv^A"^) - (v^A"^u)(uv^A~^)l 

l-v^A”^u L 


or 


I 


which is the same as the right-hand side. 
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APPENDIX B 


PROGRAM SOURCE LISTING 


MTVO UT 
PURPOSE 

IDENTIFY OUTLIERS IN MULTIVARIATE DATA 


METHOD 

LET VECTORS X( 1) THRU X( N) BE CBSERVATICNS FROM A NP-VARIATE 
NORMAL DISTRIBUTION 

THEN 

XBAR AND COVARIANCE MATRIX S ARE OBTAINED BY EQUATIONS 
XB AR = < l/N) * SUM X( J) 

S (N-1) = 5UM( X( J )-XBAR) ♦ < X ( J) -X BAR ) T RANS POS E 


WHERE 


N 
NP 
ALPH 
F ALPH 


X( J, I ) 


SAMPLE SIZE (N. LE.500 ) 

NUMBER OF VARIABLES (NP.LE.IOI 

SIGNIFICANCE LEVEL 

F-VALUE FOR ALPHA AT NP AND N-NP~L 
DEGREE S OF FREE DCM 

THE »I »TH ELEMENT CF THE »J*TH VECTOR, 
WHERE 1=1,2 ,...,NP AND J=l,2, ,N 


CALCULATE 

DELSTAR = ( ( N-1)**2«NP*FALPH) / ( ( N- NP-*l ) + NP^F AL P H ) > 


CALCULATE FOR EACH OBSERVATION VECTCR X(J) 

RR = ( X( J)-XBAR) TRANSPOSE (S)INVERSE ^ (X(J)-XBAR) 


IF RR > DELSTAR, THEN X( J) IS IDENTIFIED AS AN OUTLIER 
REFERENCE 


1. AN INTRODUCTION TO MULTIVARIATE STATISTICAL ANALYSIS, 

ANDERSON , 1965 

2. A SIMPLE TECHNIQUE FOR AUTOMATIC CCMPUTER EDITING 

OF BIODATA, NASA TN D-5275 

3. SYS/360 SCIENTIFIC SUBROUTINE PACKAGE ( 360 A-CM-0 3X ) 

PROGRAMMERS MANUAL, IBM INC., 1968 


SUBROUTINES 

FTABLE 
MPRO 
LOC 
D SINV 
DMF SD 


(IBM SSP) 

II ft 

It 11 

II 11 
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APPENDIX B 


C INPUT 

C 


C 

CARD 1 

FORMAT 

OF X-ARRAY CARDS TC BE READ IN 

(20A4) 

C 

CARD 2 




c 

COL 

1- 2 

NP (NUMBER OF VARIABLES) NP.LE.IO 

(12) 

C 

COL 

3 

BLANK 


C 

COL 

6^ 6 

ALPH (SIGNIFICANCE LEVEL) .05 OR .01 

(F3.2) 

c 

COL 

7“ 9 

BLANK 


c 

COL 

10 

VARIABLE NAME CARD INDICATOR 


c 



1 - NAME CARD FOLLCWS 


c 


BLANK - NO NAME CARD 


c 

CARD 3 

(OPTIONAL) 



C TEN FIELDS OF EIGHT CHARACTERS EACH (1CA8), WHICH MAY 

C BE USED TO ASSIGN MEANINGFUL NAMES TO THE NP VARIABLES. 

C IF COL 10 OF THE PRE VIOLS CARD IS PUNCHED, NAMES MUST 

C BE ASSIGNED FOR ALL NP VARIABLES. JDEFAULT NAMES ARE 

C ‘XI* , * X2S ... , ’X(NP) ». 

C CARDS DATA FORMATTED AS PRESCRIBED IN CARD 1 

C 

C MULTIPLE RUNS ARE PERMITTED, AS LONG AS EACH DATA DECK IS 

C PRECEDED BY APPROPRIATE CONTROL CARDS (CARDS 1, 2 AND 

C 3 ABOVE), AND IS FOLLOWED BY A CARD WITH PUNCHED 

C IN COLUMNS 1 THRU 4. 

C 

C OUTPUT 

C 1 LIST OF VECTORS WITH OUTLIERS IDENTIFIED BY ASTERISKS 

C 2 MEAN (ORIGINAL DATA) 

C 3 STANDARD DEVIATION (ORIGINAL DATA) 

C 4 MEAN OUTLIERS DELETED) 

C 5 STANDARD DEVIATION (OLTLIERS DELETED) 

C 


REALMS SUMX(IO), XBAR(IO), SD(IO), XBARl(lO), SDUIO) 

REAL=*'8 A( 55) , BllO), S(I0,10) 

REAL=«'8 RR, DELSTR, R(1C), FMT(IO), BLNK/* •/, VAR(IO) 

REALMS DEF(IO)/' XI X2 •,» X3 X4 X5 

1 X6 »,* X7 S* X8 X9 •,* XIO */ 

REAL X(501,1C), Y(50I,10) 

INTEGER SJ^R/' /f BLANK/* •/, FLAG, FLGTOT, DSWI 

INTEGER PREFMT ( 2 ) / M A4, • , *T1 , •/, CP T ( 2 ) / * YES*,* NC*/ 

INTEGER*2 IFLA G( 500) 

LOGICAL F0RMAT(8S), RPAREN/*)*/ 

EQUIVALENCE ( ^ RE FM T, F ORMA T( 1 ) ,) ,( FMT, FORMAT (9) ) ,( RPAREN ,F OP MAT (89)) 

1 DO 2 1=1,10 

VAR( I) = DEF( I ) 

2 FMT( n = BLNK 

DO 5 1=1,500 

5 IFLAG( I ) = 0 
FLGTOT = 0 
L IN E S = 1 

NSWI = 0 
DSWI = 0 
OH = 0.01 
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10 = 1 
IZ = 0 
MM = 1 
N = 1 

READ{ It 1000,E.NJD=99<3) FMT 
READ(ItIOIO) NP, ALPH, ISWI 

C TEST IF NAMES HAVE BEEN ASSIGNED TO THE VARIABLES 

IF( ISWI .GT .0 ) GO TO 10 
MM = 2 
GO TO 15 

10 REA0(1,1050) (VAR(I), 1=1, NP) 

C ^EAD IN THE DATA AS X-ARRAY 

15 READ( 1, FORMAT, END=S95) IXX, (X(N,I), I=l,NP) 

16 IF( IXX. EQ. STAR ) GO TO 19 
N = N + 1 

IF(N .GT.501) GO TO 990 
GO TO 15 

19 N = N “ 1 

C SELECT APPROPRIATE F-VALUE 

CALL FTABLE ( N ,NP , ALPH ,F A LPH ) 

C WRITE THE INPUT CONTROL INFORMATION 

WRIT E( 3, 5002 ) NP , ALPH, FALPH, OPT(MM) 

DELSTR = ((N-l)’{'^^2 ^ NP * FALPH) / (N + ( ( N-NP-1 ) + NP « FALPH)) 

20 DO 86 1= 1, 10 

SUMX( I )=0.D0 
CO 86 J= 1, 10 

86 S( I, J ) = 0 .DO 
DO 3 0 J = I , N 

C TEST FOR FLAGGED VECTORS IDENTIFIED AS OUTLIERS 

IF( IFL AG(J ) .GT .0 ) GO TO 30 
DO 25 1=1, NP 

25 SUMX( I ) = SUMX( I ) + X( J,I ) 

30 CONT INUE 

C FIND MEAN OF EACH VARIABLE 

DO 40 1=1, NP 

XBAR(I) = SUMX(I) / (N-FLGTQT) 

DO 40 J= 1,N 

Y(J, n = X( J, I ) - XBAR( I ) 

40 IF( IFLAG( J ) .GT .0 ) Y( J , I ) = C.DO 
JJ = 0 
DO 70 1=1, NP 
DO 70 K= 1, I 
DO 60 J= 1,N 

60 S(I,K) = Y(J,I) * Y(J,K) 4- S(I,K) 

C FIND STANDARD DEVIATION CF EACH VARIABLE 

SO(I) = DSQRT(S{I,I) / ( N-l-FLGTGT) ) 

JJ = JJ +1 
7C AUJ ) = S( I,K) 

IF COMPUTATIONS ARE COMPLETE, BRANCH TO 
3RINT TABLE OF MEANS AND S.D. «S 
IF(DSWI .GT .0) GO TO 20C 
CALL DSINV( A,NP,OH,IER ) 
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IF(IER) 991, 80,9S1 

C WRITE THE LIST HEADING 

30 WRITE( 3t 1015) ( VAR ( I ) , I =1 ,NP) 

DO 100 J=1»N 

CO 90 K=1,NP 
90 B(K) = Y(J,K ) 

CALL MPRD ( B,A ,R r IO,NP ,IZ f I CtNP) 

CALL MPRD (R,B ,RR, IO,NP,IZ ,IZ ,10) 

RR = RR*(N-1) 

FLAG = BLANK 

C FLAG THIS VECTOR WITH AN ASTERISK IF 

C IT IS IDENTIFIED AS AN CCTLIER 

IF(RR.GT.DELSTR) FLAG = STAR 
IF(FLAG.NE.STAR) GO TO 95 
IFLAG( J ) = 1 
FLGTDT = FLGTOT + I 

C WRITE THE DATA VECTOR AND IF IDENTIFIED AS 

C AN OUTLIER, LABEL WITH AN ASTERISK 

95 WRITE! 3, 1020) J, RR, FLAG, (MJ,K), K = 1 , NP ) 

L INES = L INES + 1 

IF(L INES .LE • 55 ) GO TO 100 

WRITE! 3, 1015) ( VAR ( I ) , 1=1, NP) 

L INES = 1 

100 CONTINUE 

WRITE! 3, 1025) ALPH 
WRITE! 3, 1030) N 
WRITE! 3, 1035) FLGTOT 
WRITE! 3, lOAO) DELSTR 

C SAVE MEANS AND S.D.'S, THEN LCCP BACK AND 

C COMPUTE NEW MEANS AND S.D. *S AFTER DELETING OUTLIERS 

DO 110 1=1, NP 

XBAR 1! I) = XBAR! I ) 
lie SDl! I) = SD( I) 

DSWI = 1 

GO TO 20 

C WRITE TABLE OF MEANS AND S.D. «S BEFORE AND 

C AFTER DELETION OF OUTLIERS 

200 WR ITE! 3, 2000 ) 

WRITE! 3, 2005) 

WRITE! 3, 2010) 

DO 210 1=1, NP 

210 WRITE! 3, 2015) VAR!I), XBARl(I), SDUI), XBAR!I), SD(I) 

IF (NSW I .EQ .1) GO TO 999 
GO TO 1 
C 

990 WRITE ( 3,5000) 

GO TO 99 99 
99 1 WRITE ! 3, 5001) 

GO TO 9999 
99 5 NSW I = 1 

GO TO 19 

999 WRITE ( 3,50C9) 



APPENDIX B 


1000 FORMAT (20A4) 

1010 FORMAT ( 12, IX, F3.2.3X, I 1)' 

1015 FORMAT <1H1,*LIST OF VECTORS WITH CUTLIERS IDENTIFIED BY ASTERISKS 
1*/ IHO, • J • , 7X, 'DELTA ' ,7X,1C(2X,A8) /) 

1020 FORMAT ( IH , 1 3 ,F 1 2 .4, 2 X , A 1 , AX , 1 OFl 0. 2) 

1025 FORMAT ( // / IHO , 1 2X , • * OUTLIER IDENTIFIED AT ‘,F4.2,' SIGNIFICANCE 
ILEVEL* ) 

1030 FORMAT ( IH ,12X, 'SAMPLE SIZE IS '13) 

1035 FORMAT ( IH ,12X, 'NO OF OUTLIERS IS ',13) 

1040 FORMAT ( IH , 12 X, 'DEL STAR = 'F10.4) 

1050 FORMAT ( lOA 8 ) 

2000 FORMAT I IH 1 , 20 X, 'M EAN AND STANDARD DEVIATION CF THE VARIABLES') 

2005 FORMAT ( IHO, 14X, 'DATA BEFORE I DE NTI F I C AT 1 0 N ' ,9 X , ' D AT A AFTER DELETI 

lON'/lH ,22X,'DF OU TL IE R S ' , 2 1 X , 'OF CUTLIERS') 

2010 FORMAT! IHO, ' VA R I AB LE S ' , 8 X , ' ME AN ' ,1 3X , ' S. D. ,' ,1 IX , ' ME AN' ,14X ,' S .D.' ) 

2015 FORMAT ( IHO, A 8 ,F 1 5 .4,F 17 . 4 ,F 1 5. 4 ,F 18 . 4 ) 

5000 FORMAT (IHl, 'SAMPLE SIZE EXCEEDS 500 - PROGRAM TERMINATED') 

5001 FORMAT (IHl, 'ERROR IN THE MATRIX INVERSION PROCESS - PROGRAM TERMI 
INATED' I 

5002 FORMAT ( IHl , '** I NP UT C ON TROL I NF ORMA TI CN '/ /IHO , ' **NUMBER OF VARIAB 
ILES IS',I4/1H ,'*-SIGNIF ICANCE LEVEL IS',F5.2/1H ,'**F-VALUE IS', 
2F16.4//1H ,'«FVARIA8LE NAME CARD; ',A4) 

5009 FORMAT (IHl, 'END OF JOB') 

9999 STOP 
END 
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SUBROUTINE FTABLE (N,M,SLtF) FTABOOOO 

Q 4 c ;}: 4 : ^: 4 : jJ: * jJ; :<c :fr 4 : 4: ^ jjc :{:* j{; ;^s s«s sjt j}: <s :i; :jc« *:{: 4s :4s :{! 4 k :Jc « jje fj A BOO 10 

=f=FTA60020 

FTABLE ^i'FTABOOBO 

=«' FT AB0040 

SUBROUTINE FTABLE SELECTS THE PPCPER VALUE CF F AT M AND N-M-L <^FTAB0050 
DEGREES OF FREEDOM FOR SIGNIFICANCE LEVELS CF 5Z OR 1^. =t^FTA80060 

♦ FTAB0070 

CALL ING PARAME TERS * FT AB0080 

4“ FT AB0090 

N = SAMPLE SIZE INTEGER «FTAB0100 

M = NUMBER OF VARIABLES M.LE.IC INTEGER *FTA60110 

SL = SIGNIFICANCE LEVEL .05 CR .01 REAL ■!<FTA80120 

F = SELECTED F-VALUE REAL «FTAB0130 

^ FT A80140 
=!=******<“ FT ABO 150 

DIMENSION TABLE( 2, 10,37) FTAB0160 

INTEGER XOF( 7) /3C, A0,6C,12G,20C,ACG,1000/, OF FTAB0170 

REAL SLT( 2 )/.0E, .01/ FTA60180 


CALL ING PARAME TERS 

N = SAMPLE SIZE 
M = NUMBER OF VARIABLE S 
SL = SIGN IF ICANCE LEVEL 
F - SELECTED F-VALUE 


M.LE. 1C 
,05 CR .01 


INTEGER 

INTEGER 

REAL 

REAL 


4t 4; 4i ^}: 4c 4: 4: 4: * 4t # 4: 4c 4t 4: * 4; 4= * 4t 4s 4; 4: 4c 4c 4c 4c 4c 4< 4t ^ 4= * 4^ 4= 4= ^ ^ ^ F T A B 0 1 5 0 

DIMENSION TABL E( 2, 10, 37) FTAB0160 

INTEGER XDFI 7) /3C, 40, 60,120, 20C, ACC, 1000/, OF FTAB0170 

REAL SLT( 2 )/.05, .01/ FTAB018C 

REAL TABL H 20 C ) / 1 6 1 . 4 , 4052 . , 1 99. 5 ,4 999. 5, 215. 7, 540 3. ,22 4.6,5625., FT ABO 190 


1230. 2t 

5764., 

234.0, 

5859. , 

236 

..8 

, 5928. , 

236 

I. 9 

,5982 . 

,240 

.5, 

60 22 . 

26056 

18.51, 

98 

.50, 

19. 

00, 

99. 

CC 

tl9 

• 16, 

99. 

17 

,19 

.25 

,99. 

25, 

19 

.30 

319.33, 

99.33, 

19 

.35, 

99. 

36, 

19. 

37 

,99 

.37, 

19. 

38 

,99 

.39 

,19. 

40, 

99 

.40 

434.12, 

9.55, 

30 

.82, 

9. 

23, 

29. 

46 

, 9 

.12 , 

28. 

71 . 

, 9 

.0 1 

,28. 

24, 

8 

.94 

5 

8 .89, 

27.67, 

8 

.85, 

27. 

49, 

8. 

81 

,27 

.35, 

8. 

79 

,27 

C\J 

, 7. 

71, 

21 

.20 

618.00, 

6.59, 

16 

.69, 

6. 

39, 

15. 

98 

, 6 

.26, 

15. 

52 . 

r 6 

.16 

,15. 

21, 

6 

.09 

7 

6 .04, 

14.80, 

6 

.00, 

14. 

66, 

5. 

96 

,14 

.55, 

6. 

61 - 

rl6 

.26 

, 5. 

79, 

13 

.27 

812 .06, 

5. 19, 

11 

.39, 

5. 

05, 

10. 

97 

, 4 

.95, 

10. 

67 

, 4 

.38 

, 10 . 

46 , 

4 

.82 

9 

4.77, 

10.16, 

4 

.74, 

10. 

05, 

5. 

99 

,13 

.75, 

5. 

14 : 

riO 

.92 

, 4 • 

76, 

9 

.78 

1 

9.15, 

4.39, 

8 

.75, 

4. 

28, 

8. 

47 

, 4 

.21 , 

8. 

26, 

r 4 

.15 

, 8 • 

10, 

4 

. 10 

2 

4.06, 

7.87, 

5 

.59, 

12. 

25, 

4. 

74 

, 9 

.55, 

4. 

35, 

) 8 

.45 

, 4 . 

12, 

7 

.85 

3 

7.46, 

3.87, 

7 

.19, 

3. 

79, 

6. 

99 

, 3 

. 73 , 

6. 

84, 

) 3 

.68 

, 6 • 

72, 

3 

.64 

4 

5.32, 

11.26, 

4 

.46, 

8. 

65, 

4. 

07 

, 7 

. 59, 

3. 

84 , 

r 7 

.01 

, 3. 

59 , 

6 

.63 

5 

6.37, 

3.50, 

6 

.18, 

3. 

44, 

6. 

C3 

, 3 

.39, 

5. 

91 , 

) 3 

.55 

t 5. 

3 1, 

5 

.12 

6 

4.26, 

8.02, 

3 

.86, 

6. 

99, 

3. 

63 

, 6 

.42 , 

3. 

48, 

r 6 

.06 

, 3. 

37, 

5 

.SO 

7 

5.6 1, 

3.23, 

5 

.47, 

3. 

18, 

5. 

35 

, 3 

.14, 

5. 

26, 

) 4 

.96 

, 10. 

04, 

4 

.10 

8 

3.71, 

6.55, 

3 

.48, 

5. 

99, 

3. 

33 

, 5 

.64 r 

3. 

22, 

r 5 

.39 

, 3. 

14, 

5 

.20 

9 

5.06, 

3.02, 

4 

.94, 

2. 

98, 

4. 

85/ 










REAL 

TABL2( 

200) 

/ 4 

. 84 

,9. 

65 

,3. 

98 ,7 

.21 

,3. 

59 

,6 . 

22 ,3 

.36 

,5 

.67 

15 

.32,3 

.09, 5. 

07 

,3.01,4 

. 89 

,2. 

95 

,4. 

74,2 

. 90 

,4. 

63 

,2 . 

85 ,4 

.54 

,4 

.75 

23 

.89,6 

.93, 3. 

49 

,5.9 

5, 3 

. 26 

,5. 

41 

,3. 

11 ,5 

. 06 

,3. 

00 

,4 . 

8 2,2 

.9 1 

,4 

.64 

34 

.50, 2 

. 80, 4. 

39 

,2.75,4 

. 30 

, 4 • 

67 

,9. 

07 ,3 

.81 

,6. 

70 

,3 . 

41,3 

.74 

,3 

. 18 

43 

.03,4 

. 86, 2. 

92 

,4.62,2 

. 83 

,4. 

44 

, 2 . 

77 ,4 

.30 

,2. 

71 

,4 . 

19 ,2 

.67 


.10^ 

58 

.86, 3 

.74, 6. 

51 

,3.34,5 

. 56 

,3. 

11 

,5. 

04 ,2 

. 96 

,4 • 

69 

,2 . 

85,4 

• 46 

,2 

.76 

62 

.70, 4 

. 14, 2. 

65 

,4.0 

3,2 

. 60 

,3. 

94 

, 4 • 

54,8 

.68 

,3. 

68 

,6 . 

36,3 

.29 

,5 

.42. 

74 

.89, 2 

.90, 4. 

56 

,2.79,4 

. 32 

,2. 

71 

, 4 . 

14 ,2 

. 64 

,4. 

00 

,2 . 

39 ,3 

.89 

,2 

.54: 


FT A60370 


FT4B0420 


8 4.49, 8.5 3, 3.63,6.23,3.24,5.29,3.01 ,4.77 ,2. 85 ,4 .44,2.74,4.20,2.66, FT AB0460 
94 .03, 2 . 59, 3. 89,2 . 54,3. 78 ,2.49 ,3. 69,4.45 ,8. 40,3 .59 ,6 . 1 1 ,3 .20 , 5 . 18, FT AB0470 
12.96, 4.67, 2.81 ,4.34, 2. 70,4. 1C,2. 61 ,3 .93 ,2. 55 ,3 .79 ,2 .49 ,3 .68 , 2 .45, FT AB0480 
23. 59, 4. 41, 8. 29, 3. 55, 6. 01, 3. 16, 5. 09, 2. 93 , 4. 58, 2. 77, 4. 25, 2. 66, 4. 01, FT A3 0490 
32. 58, 3. 84, 2. 51, 3. 71, 2. 46, 3. 60,2.41 ,3 . 51 , 4. 38, 8. 18, 3. 52, 5 . 9 3, 3. 13, FT A30500 
45.01, 2 .90, 4. 50,2.74, 4. 17,2. 63,3. 94,2 . 54,3. 77 ,2 .48 ,3 .63 ,2 .42, 3 .52, FT A80510 
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52.38,3«A3f4.35T8.l0»3*49t5.85t3.10f4.94i2.37T^» 
6 3. 87t 2 ,5 If 3,70, 2 .45,3* 56 f2. 3 9, 3. 46 ,2 .35 ,3 • 37/ 

PEAL TABL3(200>/ 4 . 32 t 8 . C2 , 3 . 47 »5 . 78 1 3 . 07 ,4 . 

14.04, 2. 57f 3. 31 f2 .49,3.64t2. 42,3. 5L ,2 .37 ,3.40t2 . 

23-44, 5 .72, 3.05,4. 82, 2. 82,4.31 ,2. 66,3 .99 ,2, 55 ,3 . 

3 3.45, 2.34, 3.35,2.3 0,3. 26,4. 2 8,7. 83 ,3 .42 ,5. 66,3 . 

42.64, 3.94, 2.53, 3.71, 2. 44,3. 54,2. 37, 3 .41 ,2. 32, 3 . 
57.82,3. 40, 5.61,3.01,4.72,2. 78,4. 22 ,2 . 62 , 3. 90 ,2. 

62.36.3.36, 2.30,3.26,2.25,3.17,4.24,7.77 ,3.39,5 . 

74. 18, 2. 60, 3. 85, 2. 49, 3. 63, 2.40,3.46 ,2 . 34 ,3. 32, 2. 
3 4.23,7.7 2, 3.37,5.5 3,2.98,4.64,2. 74,4.14,2. 59,3 . 

9 3 .42, 2.32, 3. 29, 2. 27, 3. 18, 2. 22, 3. C9, 4. 21, 7. 68, 3. 

12.7 3, 4.11, 2. 57, 3. 7 8, 2. 46, 3. 56, 2. 37 , 3 . 39, 2. 31 , 3 . 

23.06,4.20, 7.64,3.34,5.45,2.95,4.57 ,2.71 ,4. 07 ,2 • 

32.36, 3.36, 2. 29, 3. 23, 2. 24, 3. 12, 2. 19, 3. 03 ,4. 18, 7. 
44.54, 2.70, 4.04,2.55,3.73,2.43,3. 50, 2. 35, 3. 33, 2 . 

5 2.18, 3.00, 4. 17,7. 5 6,3.32,5.39,2. 92 ,4.51 ,2. 59,4. 

6 3.47,2.3 3, 3. 30, 2. 27, 3. 17, 2. 21 , 3. 07, 2. 16, 2. 98/ 

REAL TABL4( 140)/ 4.08,7. 31,3.23,5.18,2. 84,4. 

13. 51, 2. 3 4, 3. 29, 2. 25, 3. 12, 2. 18, 2. 99, 2 . 12, 2. 89, 2 . 

23. 15, 4.9 8, 2.76,4.13,2. 53, 3. 65, 2. 37, 3 .34 ,2. 25,3 • 

3 2.82, 2.04, 2. 72 ,1.99, 2. 63, 3. 92, 6. 85, 3. 07, 4. 79, 2. 
42. 29, 3 . 17, 2. 17, 2 .96, 2. 09,2. 79,2. 02 ,2 .66 ,1. 96 ,2 . 

56.76, 3.0 4, 4. 71, 2. 65, 3. 88, 2. 41, 3. 41 ,2 .26, 3. 11,2 . 

61.98.2.6 0, 1.92,2.50,1.87 ,2.41,3. 86 ,6.70,3. 02,4. 

73.36, 2 . 23, 3.06,2.12, 2. 85,2. 03,2.69, 1 .96, 2. 55,1 

83.85.6.6 6, 3. 00, 4. 62, 2. 61, 3. 80,2.38, 3 .34 , 2. 22,3 

9 2.66, 1 .9 5, 2. 53, I .89, 2.43, 1. 84,2. 34 ,3 .84 ,6. 63 ,3 

12 .37, 3 .32, 2. 21 ,3 .02, 2 . 10,2 . 8C,2. 01 ,2 . 64 ,1. 94 ,2 

22.32/ 

EQUIVALENCE ( T AB LE ( 1 , 1 , 1 ) , TABL 1 ( 1) ) , (TABLEd, 
EQUIVALENCE ( T AB LE ( L , 1 , 2 1) , TAB L 3 ( 1 ) ) , (TABLE(1, 
IF(M .GT .10) G3 TO 590 


43,2.71,4.10,2.60, 


87.2 

32.3 
76 ,2 
,03 ,4 

30.2 
51 ,3 

57.2 
28 ,3 
82 ,2 
35 ,5 

26.2 
56 ,3 
60 ,3 
28 ,3 
02 ,2 


.84,4 
.31,4 
.46 ,3 
.76 ,2 
.27 ,3 
.67 ,2 
.99 ,4 
.22 ,2 
.47 ,3 
.49 ,2 
.25 ,3 
.75 ,2 
.33 ,5 
.20 ,2 
.53,3 


.37,2 
.30,7 
.59,2 
.80,4 
.21,4 
.42,3 
. 68, 2 
.24,3 
.5<^,2 
.96,4 
.15,2 
45,3 

42.2 

22.3 
70,2 


.50, 
.76, 
.13, 
.39, 
.60, 
. 20 , 


.09 , 
.42i 


31.2 

08 ,2 

12 .2 

68,3 

56 .1 

14.2 

66.2 
90 ,2 
04,2 
00 ,4 
51,1 


,61 ,3 
.80 ,4 
,17 ,2 
,95,2 
.91,2 
.90,2 
.62 ,3 
.46 ,1 
. 10 ,2 
.6 1,2 
.83 ,2 


.83,2 
.00,7 
.95, 2 
.45,3 
.47,3 
.05, 2 
. 83,2 
.85, 2 
.82,2 
.60 , 3 
41,1 


1,11) ,TABL2ll) 
1 ,31 ) ,T ABL4 ( 1 ) 


DF = N-M- 1 

IF(SL .EQ .SLT( 1 ) ) GO TO 10 
IF( SL .EG .SLT< 2 ) ) GO TO 1 5 
GO TO 591 
10 L = 1 

GO TO 16 

15 L = 2 

16 IF(DF.LE.30) GO TO 40 

DO 20 1=2,7 

1F(DF~XDF( II) 50, 3 C, 20 

20 CONTINUE 

F = TABLEIL ,M, 36) + 

RETURN 

30 DF = 1+29 

40 F = TABLEtL ,M, DF ) 

50 F = TABL E( L ,M, 1+29) + ((l./OF - l./XDF{I)) / {I./XDF(I-l) ~ 1 
1(1))) - ( TABLE (L ,M , 1 + 28) - TAB LE I L , M , I +29) ) 

RETURN 


( ( 1. /N) /.OCl ) * (TABLE (L,M,36)-T ABLE(L,M, 


PT AB0520 
FT A30530 
.68, FTAB0540 
.95, FTA30550 

.40, FTA805&0 

.26, FTAB0570 
.26, FTAB0580 
FT AB0590 
FT AB0600 
FT AB0610 
FT AB0620 
FTAB0630 
FT A0O64O 
.53, FTAB0650 
.93, FTA0O66O 
FT AB0670 
FT AB0680 
FT AB0690 
FT AB0700 
FTAB07L0 
FT AB0720 
FT AB0730 
FT AB0740 
FT AB0750 
FT AB0760 
FT AB0770 
FT AB0780 
.78, FTAB0790 
.83, FTAB0800 
FT AB08 10 
) FTAB0820 

) FTA80830 

FT A60840 
FT AB0850 
FT AB0660 
FT AB0870 
FT AB0880 
FT AB0890 
FT A30900 
FT ABC9 10 
FT AB09 20 
FT AB09 30 
FT AB0940 
FT AB0950 
37 ) ) FT AB0960 
FT AB0970 
FT AB0980 
FT AB0990 
FT A3 1000 
./XDFFT ABlOlO 
FT A31020 
FT AB 10 30 


.45, 

.08, 

. 10 , 

.48, 

.89, 

.73, 

.39, 

.37, 

. 02 , 


590 WRITE! 3,990) 

F = 2.32 

RETURN 

591 WRITE(3,991) SL 
SL = 0.01 
GO TO 15 

990 FORMAT ( IH 1 , *N UMBER OF VARIABLES > 10. 

991 ^F 0 RMAT^t lHld*^*SIGNI^ LEVEL •,F6.3,« IS NOT ACCEPTABLE. 

lEL IS SET TO .Cl . * ) 

END 


F HAS BEEN S ET T 0 A 


FT A81040 
FT AB 10 50 
FT AB 1060 
FT AB1070 
FT AB1080 
FT AB 1090 
DUMMFT AB 1100 
FT ABlllO 
LEVFTAB1120 
FT A31130 
FT AB1140 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE MPRO 
PURPOSE 

MULTIPLY TWO MATRICES TO FORM A RESULTANT MATRIX 
USAGE 

CALL MPRD?AtB,R,N»M,MSA,MSBTL< 


DESCRIP TION 


A - 
B - 
R - 
N - 
M - 
MSA 


MSB 
L - 


NAME 

NAME 

NAME 


OF PARAMETERS 
OF FIRST INPLT MATRIX 
OF SECOND INPLT MATRI X 
OF OUTPUT MA TRI X 


NUMBER OF ROWS IN A AND R 

NUMBER OF COLUMNS IN A AND ROWS IN B 

- ONE DIGIT NUMBER FOR STORAGE MODE OF 


MATRIX A 


0 - 
1 - 
2 - 
- SAME 
NUMBER 


GENERAL 
SVMME TRIC 
D lAGONAL 
AS MSA EXCEPT 
OF COLUMNS IN 


FOR MATRIX 
B AND R 


REMARKS 

MATRIX R CANNOT BE IN THE SAME LCCATIGN AS MATRICES A OR 
NUMBER OF COLUMNS OF MATRIX A MUST BE EQUAL TO NUMBER OF 
OF MATRI X B 


SUBROUTINES 

LOC 


AND FUNCTION SUBPROGRAMS REQUIRED 


METHOD 

THE M BY L MATRIX B IS PRE 
AND THE RESULT IS STORED I 
ROW INTO COLUMN PRODUCT, 
THE FOLLOWING TABLE SHOWS 
MATRIX FOR ALL COMBINATION 
A 

GENERAL 
GENERAL 
GENERAL 
SYMME TRIC 
SVMME TRIC 
SYMME TRIC 
D lAGONAL 
DIAGONAL 
D lAGONAL 


MULTIPLIED 
N THE N BY 


BY THE N 
L MATRIX 


BY 

R. 


M MATRIX 
THIS IS 


THE STORAGE MODE OF 
MATRICES 


THE OUTPUT 


S OF 


I NPUT 
B 

GE NERAL 
SYMMETRI C 
DI AGONAL 
GE NEPAL 
SYMMETRIC 
DI AGONAL 
GENERAL 
SYMMETRI C 
DI AGONAL 


GENERAL 

GENERAL 

GENERAL 

GENERAL 

GENERAL 

GENERAL 

GENERAL 

GENERAL 

DIAGONAL 


SUBROUT INE M PR DC A , B ,R , N , M , M SA , MSB , L) 


MPRD 
, . ..MPRO 
MPRO 
MPRD 
MPRO 
MPRO 
MPRO 
MPRO 
MPRD 
MPRD 
MPRD 
MPRD 
MPRD 
MPRD 
MPRD 
MPRD 
MPRD 
MPRD 
MPRD 
MPRD 
MPRO 
MPRD 
MPRD 
MPRO 
MPRD 
B MPRD 
ROWMPRD 
MPRD 
MPRD 
MPRD 
MPRO 
MPRD 
MPRD 
A MPRO 
A MPRD 
MPRD 
MPRD 
MPRD 
MPRD 
MPRO 
MPRD 
MPRO 
MPRO 
MPRD 
MPRD 
MPRD 
MPRD 
MPRD 
MPRD 
. . .MPRD 
MPRO 
MPRO 


001 

002 

003 

004 

005 

006 

007 

008 

009 

010 

011 

012 

013 

014 

015 

016 

017 

018 
019 
0 20 
021 
022 

023 

024 

025 

026 

027 

028 

029 

030 

031 

032 

033 

034 

035 

036 

037 

038 

039 

040 

041 

042 

043 

044 

045 

046 

047 

048 

049 

050 

051 

052 
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DOUBLE PRECISION A, B, R 



DIMENSION At 1) ,B( 1),R( 1) 

MPRD 

053 



MPRD 

054 


SPECIAL CASE FOR DIAGONAL BY DIAGONAL 

MPRD 

055 



MPRD 

056 


MS=MSA* 10+MSB 

MPRD 

057 


IF(MS-22) 30, 1C, 30 

MPRD 

058 

10 

DO 20 1=1, N 

MPRD 

059 

20 

R( I ) = A( I )*B( n 

MPRD 

060 


RETURN 

MPRD 

061 



MPRD 

062 


ALL OTHER CASES 

MPRD 

063 



MPRD 

064 

30 

IR=1 

MPRD 

065 


DO 90 K=1,L 

MPRD 

066 


DO 90 J= 1,N 

MPRD 

067 


R( IR )=0 

MPRD 

068 


DO 80 1=1, M 

MPRD 

069 


IF(MS) 40, 60, 4C 

MPRD 

070 

40 

CALL LOCtJ, I,IA,N,M,MSA) 

MPRD 

071 


CALL LOC( I,K,IB,M,L,MSB> 

MPRD 

072 


IF(IA) 50,8C,5C 

MPRD 

073 

50 

IF(IB) 70,80,7C 

MPRD 

074 

60 

IA=N’«=( I- 1) +J 

MPRD 

075 


IB=M4^{ K- 1) +I 

MPRD 

076 

70 

R ( IR )-R ( IR )+A( lA )*B( IB ) 

MPRD 

077 

30 

CONT INUE 

MPRD 

078 

90 

IR= IR + 1 

MPRD 

079 


RETURN 

MPRD 

080 


END 

MPRD 

081 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 





LOC 
• LOC 

001 

002 




LOC 

003 

SUBROUTINE LOC 


LOC 

004 




LOC 

005 

PURPOSE 



LOC 

006 

COMPUTE A VECTOR SUBSCRIPT FOR AN ELEMENT IN A MATRIX OF 


LOC 

007 

SPECIFIED STORAGE MODE 


LOC 

008 




LOC 

009 

USAGE 



LOC 

010 

CALL 

LOC 2It J, IR tNtM ,MS< 


LOC 

Oil 




LOC 

012 

DESCRIPT 

ION OF PARAMETERS 


LOC 

013 

I 

ROH NUMBER OF ELEMENT 


LOC 

014 

J 

COLUMN NUMBER OF ELEMENT 


LOC 

015 

IR - 

RESULTANT VECTOR SUBSCRIPT 


LOC 

016 

N 

NUMBER OF ROkS IN MATRIX 


LOC 

017 

M 

NUMBER OF COLUMNS IN MATRIX 


LOC 

018 

MS - 

ONE DIGIT NUMBER FOR STORAGE MCDE OF MATRIX^ 


LOC 

019 


0 - GENERAL 


LOC 

020 


1 - S VMME TRIG 


LOC 

021 


2 - DIAGONAL 


LOC 

022 




LOC 

023 

REM ARKS 



LOC 

024 

NONE 



LOC 

025 




LOC 

026 

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 


LOC 

027 

NONE 



LOC 

028 




LOC 

029 

M ETHOD 



LOC 

0 30 

MS#0 

SUBSCRIPT IS COMPUTED F CR A MATRIX WITH N4^M ELEMENTS 

LOC 

031 


IN STORAGE ^GENERAL MATRIX< 


LOC 

032 

MS#1 

SUBSCRIPT IS COMPUTED FOR A MATRIX WITH N=4'^NCl</2 

IN 

LOC 

033 


STORAGE SUPPER TRIANGLE GF SYMMETRIC MATRIX<. IF 


LOC 

034 


ELEMENT IS IN LOWER TRIANGULAR PORTION, SUBSCRIPT 

IS 

LOC 

035 


CORRESPONDING ELEMENT IN UPPER TRIANGLE. 


LOC 

036 

MS#2 

SUBSCRIPT IS COMPUTED FOR A MATRIX WITH N ELEMENTS 


LOC 

037 


IN STORAGE IDI AGONAL ELEMENTS OF DIAGONAL MATRIX<. 


LOC 

038 


IF ELEMENT IS NCT ON DIAGONAL ?ANO THEREFORE NOT IN 

LOC 

039 


STGRAGE<, IR IS SET TC ZERC. 


LOC 

040 


LOC 041 




LOC 

043 


SUBROUTINE LOC ( I , J , I R , N , M , M S ) 

LOC 

044 



LOC 

045 


IX= I 

LOC 

046 


JX=J 

LOC 

047 


IF(MS-l) 10,20,3C 

LOC 

048 

10 

IRX=N*( JX- l)+i X 

LOC 

049 


GO TO 36 

LOC 

050 

20 

IF( IX-JX) 22,24, 24 

LOC 

051 

22 

IRX= IX + { JX^JX- JX)/2 

LOC 

052 


GO TO 36 

LOC 

053 

24 

IRX=JX + ( IX* IX- IX)/2 

LOC 

054 


GO TO 36 

LOC 

055 

30 

IRX=0 

LOC 

056 


IF( IX-JXJ 36,32, 36 

LOC 

0 57 

32 

IRX= IX 

LOC 

058 

36 

IR= IRX 

LOC 

059 


RETURN 

LOC 

060 


END 

LOC 

061 
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APPENDIX B 


c 

c 


SINV 

SINV 

SINV 

SUBROUTINE DSINV 

S INV 

PURPOSE SINV 

INVERT A GIVEN SYMMETRIC POSITIVE DEFINITE MATRIX SINV 

S INV 

USAGE SINV 

CALL DSINVIA ,NtEPS,IER) 

SINV 

DESCRIPTION OF PARAME TE R S SINV 

A - DOUBLE PRECISION UPPER TRIANGULAR PART OF GIVEN 

SYMMETRIC POSITIVE DEFINITE N BY N COEFFICIENT 
MATRI X. 

ON RETURN A CONTAINS THE RESULTANT UPPER SINV 

TRIANGULAR MATRIX IN DOUBLE PRECISION. 

N - THE NUMBER OF ROWS (COLUMNS) IN GIVEN MATRIX. SINV 

EPS - SINGLE PRECISION INPUT CONSTANT WHICH IS USED 

AS RELATIVE TOLERANCE FCR TE,ST ON LOSS OF 
SIGNIFICANCE . 

lER - RESULTING ERROR PARAMETER C CDED AS FOLLOWS SINV 

IER = 0 - NO ERROR S INV 

IER=-1 - NO RESULT BECAUSE CF WRONG INPUT PARAME- SINV 
TER N OR BECAUSE S CME RADICAND IS NON- SINV 

POSITIVE (MATRIX A IS NOT POSITIVE SINV 

DEFINITE, POSSIBLY DUE TC LOSS OF SIGNI- SINV 
FICANCE) SINV 

IER=K - WARNING WHICH INDICATES LOSS CF SIGNIFI- SINV 
CANCE. THE RADICAND FORMED AT FACTQRIZA- SINV 
TION STEP K+1 WAS STILL POSITIVE BUT NO SINV 
LONGER GREATER THAN AB S ( E PS ^ A ( K+ 1 , K+ 1 )J . SINV 

SINV 

REMARKS SINV 

THE UPPER TRIANGULAR PART OF GIVEN MATRIX IS ASSUMED TO BE SINV 
STORED COLUMNWISE IN N*(N+U/2 SUCCESSIVE STORAGE LOC AT I ONS . S I NV 
IN THE SAME STORAGE LOCATIONS THE RESULTING UPPER TRIANGU- SINV 
LAR MATRIX IS STORED COLUMNWISE TCC. SINV 

THE PROCEDURE GIVES RESULTS IF N IS GREATER THAN 0 AND ALL SINV 
CALCULATED RADICANDS ARE POSITIVE. SINV 

S INV 

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED SINV 

DMFSO 

S INV 

METHOD SINV 


SOLUTIDM IS DONE USING FACTORIZATION BY SUBROUTINE DMFSD. 

S INV 

SINV 

S INV 

SUBROUTINE D SI N V( A ,N ,E P S , I E R ) 

S INV 
SINV 


10 

20 

30 

50 

60 

70 

80 

90 

110 

120 


150 

170 


200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

420 
4 30 

450 

460 

470 

490 

500 
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APPENDIX B 


DIMENSION A( 1) 

DOUBLE PRECISION A, DIN, WORK 
C 

C FACTORIZE GIVEN MATRIX BY MEANS OF SUBROUTINE DMFSD 

C A = TRANSP3 SE ( T ) * T 

CALL DMFSDt A,\J ,EPS,IER ) 

IF( lER 5 9, 1, 1 
C 

C INVERT UPPER TRIANGULAR MATRIX T 

C PREPARE INVER SI ON- LOOP 

1 IPIV=N*(N+l)/2 
IND= IP IV 

C 

C INIT lAL IZE INVERSION-LOOP 

DO 6 1=1, N 
DIN= l.DO/A( IPI V) 

A( IP IV )=DIN 
MIN=N 
KEND=I-1 
LANF=N-K END 
IF(KEND) 5, 5,2 

2 J= IND 
C 

C INIT lAL IZE ROW-LOOP 

DO 4 K= ItKEND 
WGRK = 0 .DO 
MIN=MIN- 1 
LHOR=IP IV 
LVER=J 
C 

C START INNER LOOP 

DO 3 L = LANF,MIN 
LVER=LVER+l 
LHOR = LHOR+L 

3 WORK=WORK + A( L VER )=^=A( LHOR ) 

C END OF INNER LOOP 

C 

A( J J=-WORK^D IN 

4 J = J-MIN 

C END OF ROW- LOOP 

C 

5 IP IV= IP IV-M IN 

6 IND=IND-l 

C END OF INVERSION-LOOP 

C 

C CALCULATE INVERSE(A) BY MEANS OF INVERSE(T) 

C INVERSE(A) = INVERSE(T) ♦ TRA N SPO SE ( I NVE RSE (T ) ) 

C INITIALIZE MULTIPLICATION-LOOP 

DO 8 1=1, N 
IPIV=IP IV+I 
J= IP IV 
C 


SINV 510 

SINV 530 

SINV 550 

SINV 570 
SINV 580 
SINV 590 
SINV 600 
SINV 610 
SINV 620 
SINV 630 
SINV 640 
SINV 650 

SINV 670 
SINV 680 
SINV 690 
SINV 700 
SINV 710 
SINV 720 
SINV 730 
SINV 740 
SINV 750 
SINV 760 
SINV 770 
SINV 780 
SINV 790 
SINV 800 
SINV 810 
SINV 820 
SINV 830 
SINV 840 

SINV 860 
SINV 870 
SINV 880 
SINV 890 
SINV 900 
SINV 910 
SINV 920 
SINV 930 
SINV 940 
SINV 950 
SINV 960 
SINV 970 
SINV 980 
SINV 990 
SINV 1000 
SINV 1010 
SINV 1020 
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APPENDIX B 


C INITIAL IZE ROW-LOOP 

DO 8 K=I,N 

WORK=O.DO 
LHOR=J 
C 

C START INNER LOOP 

DO 7 L = K,N 

LVER = L HOR+K- I 
WORK = WORK+A(LHOR )*A( LVER) 

7 LHOR = LHOR+L 

C END OF INNER LOOP 

C 

A(J )=WORK 

8 J = J+K 

C END OF ROW- AND MUL TI PL IC ATI CN-LOCP 

C 

9 RETURN 
END 


S INV 1030 
SINV lOAO 
SINV 10 50 
S INV 1060 
SINV 1070 
SINV 1080 
SINV 1090 
SINV 1100 

SINV 1120 
SINV 1130 
SINVllAO 
SINV 1150 
S INV 1160 
SINV 1170 
S INV 1180 
SINV 1190 
SINV 1200 


I 
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C MFSO 10 

C MFSO 20 

C MFSD 30 

C SUBROUTINE DMFSD 

C MFSD 50 

C PURPOSE MFSD 60 

C FACTOR ^ GIVEN SYMMETRIC POSITIVE DEFINITE MATRIX MFSD 70 

C MFSD 80 

C USAGE MFSD 90 

C CALL DMFSDI A,N,EPS,IER) 

C MFSD 110 

C DESCRIPTION OF PARAMETERS MFSD 120 

C A - DOUBLE PRECISION UPPER TRIANGULAR PART OF GIVEN 

C SYMMETRIC POSITIVE DEFINITE N BY N COEFFICIENT 

C MATRIX. 

C ON RETURN A CONTAINS THE RESULTANT UPPER MFSD 150 

C TRIANGULAR MATRIX IN DCLBLE PRECISION. 

C N - THE NUMBER OF ROV^S (CCLUMNSI IN GIVEN MATRIX. MFSD 170 

C EPS - SINGLE PRECISION INPUT CONSTANT WHICH IS USED 

C AS RELATIVE TCLERANCE FOR T^ST ON LOSS OF 

C SIGNIFICANCE . 

C lER - RESULTING ERROR PARAMETER C COED AS FOLLOWS MFSD 200 

C IER = 0 - NO ERROR MFSO 210 

C IER=-1 - NO RESULT BECAUSE CF WRONG INPUT PARAME- MFSD 220 

C TER N OR BECAUSE S CME RADICAND IS NON- MFSO 230 

C POSITIVE (MATRIX A IS NOT POSITIVE MFSD 2A0 

C DEFINITE, POSSIBLY DUE TC LOSS OF SIGNI- MFSO 250 

C FICANCE) MFSO 260 

C IER=K - WARNING WHICH INDICATES LOSS OF SIGNIFI- MFSD 270 

C CANCE. THE RADICAND FORMED AT FACTORIZA- MFSD 280 

C TION STEP K+1 WAS STILL POSITIVE BUT NO MFSD 290 

C LONGER GREATER THAN AB S ( E PS * A ( K+ 1 , K+ 1) ) . MFSO 300 

C MFSD 310 

C REMARKS MFSD 320 

C THE UPPER TRIANGULAR PART OF GIVEN MATRIX IS ASSUMED TO BE MFSD 330 

C STORED COLUMNWISE IN N*(N+l)/2 SUCCESSIVE STORAGE LOCAT I GNS . MFSO 340 

C IN THE SAME STORAGE LOCATIONS THE RESULTING UPPER TRIANGU- MFSO 350 

C LAR MATRIX IS STORED COLUMNWISE TOC, MFSD 360 

C THE PROCEDURE GIVES RESULTS IF N IS GREATER THAN 0 AND ALL MFSD 370 

C CALCULATED RADICANOS ARE POSITIVE. MFSD 380 

C THE PRODUCT OF RETURNED DIAGONAL TERMS IS EQUAL TO THE MFSD 390 

C SQUARE-ROOT OF THE DETERMINANT OF THE GIVEN MATRIX. MFSD 400 

C MFSD 410 

C SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED MFSD 420 

C NONE MFSD 430 

C MFSD 440 

C METHOD MFSD 450 

C SOLUTIOVI IS DONE USING THE SQUARE-ROCT METHOD OF CHOLESKY. MFSD 460 

C THE GIVEN MATRIX IS REPRESENTED AS PRODUCT OF TWO T R I ANGUL ARMFSD 470 

C MATRICES, WHERE THE LEFT HAND FACTOR IS THE TRANSPOSE OF MFSD 480 

C THE RETURNED RIGHT HAND FACTOR. MFSD 490 

C MFSD 510 
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C 

c 

SUBROUTINE DMF SD ( A , N tE P S , I E R ) 


DIMENSION A( IJ 

DOUBLE PRECISION DPIV,OSUM,A 

TEST ON WRONG INPUT PARAMETER N 
IF{N-n 12, 1,1 

1 IER=0 

INITIALIZE DIAGONAL-LOOP 
KP IV = 0 
DO 11 K=1,N 
KP IV = KP IV+K 
1ND=KP IV 
LEND=K-1 

CALCULATE TOLERANCE 
TOL= AB S ( EP S* SNl GL ( A ( KP I V ) ) ) 

START FACTORIZATION-LOOP OVER K 
DO 11 I=K,N 
CSUM=0 .DO 
IF(L END) 2, A, 2 

START INNER LOOP 

2 DO 3 L= 1,LEND 
LANF=KP IV-L 

L IND= IND-L 

3 DSUM=DSUM+A(LANF j*A( LIND) 

END OF INNER LOOP 

TRANSFORM ELEMENT A ( I NO ) 

A DSUM=A( IND)-DSUM 
IF( I-K ) 10, 5, 1C 

TEST FOR NEGATIVE PIVOT ELEMENT 

5 IF( SNGL { DSUM )- TOL ) 6,6,9 

6 IF(DSUM) 12,12,7 

7 IF( lER ) 8, 8, 9 

8 IER=K-1 

COMPUTE P IVOT ELEMENT 

9 DP IV = DSQRT( DSUM) 

A(KP IV )=DP IV 

DP IV= l.DO/DP IV 
GO TO 11 

CALCULATE TERMS IN ROW 

10 A( IND) = DSUM*DP IV 

1 1 IND= IND+I 

END OF DIAGONAL-LOOP 
RETURN 

12 IER=-1 
RETURN 
END 


MFSD 500 
MFSD 520 

MFSD 5A0 
MFSD 550 
MFSD 560 

MFSD 580 
MFSD 590 
MFSD 600 
MFSD 610 
MFSD 620 
MFSD 630 
MFSD 6A0 
MFSD 650 
MFSD 660 
MFSD 670 
MFSD 680 
MFSD 690 
MFSD 700 

MFSD 720 

TH RCW MFSD 730 

MFSD 740 
MFSD 750 
MFSD 760 
MFSD 770 
MFSD 780 
MFSD 790 
MFSD 800 
MFSD 810 

MFSD 830 
MFSD 8A0 
MFSD 850 

MFSD 870 
MFSD 880 

AND FCR LOSS CF SIGNIFICANCE MFSD 890 

MFSD 900 
MFSD 910 
MFSD 920 
MFSD 930 
MFSD 9A0 
MFSD 950 
MFSO 960 
MFSD 970 
MFSD 980 
MFSD 990 
MFSDIOOO 
MFSDIOIO 
MFS01020 
MFSD1030 
MFSD lOAO 
MFSD1050 
MFSD1060 
MFSD1070 
MFSDlOaO 
MFS01090 
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OUTPUT FOR THE EXAMPLE 


LIST CF VECTORS WITH OUTLIERS IDENTIFIER BY ASTERISKS 


J 

DELTA 

H/R 

J 

DELTA 

H/R 

1 

0.6901 

91.12 

34 

1.13 16 

89 .42 

2 

0.02^0 

96.05 

35 

0.8436 

90.48 

3 

0.01 39 

96.32 

36 

2.7647 

8 5.05 

4 

0.492^ 

1C2.30 

37 

0.0249 

96.03 

5 

0.3994 

92.57 

38 

0.0203 

98.22 

6 

0.0446 

95.64 

39 

0.0452 

95.63 

7 

0.6192 

91 .44 

40 

1 . 58 3 1 

106 .36 

8 

0. 30 98 

93.12 

41 

0.0170 

96 .2 3 

9 

0 . 1806 

94.08 

42 

0.7866 

9C.71 

10 

0. 52 39 

9 1.90 

43 

0. 1973 

9 3.94 

11 

0.0063 

96.60 

44 

0. 2601 

93.46 

12 

0.01 73 

98.14 

46 

0. 2643 

93.43 

13 

0.3266 

101.35 

4b 

0.0440 

95.65 

14 

0. 1055 

99.55 

4 7 

0.0087 

96.50 

15 

G.0390 

95 . 74 

48 

0.4256 

lOi .94 

15 

0,0244 

98.32 

49 

2. 19 12 

107.98 

17 

0 . 14 20 

99.93 

50 

0.0963 

99.45 

18 

0.68 30 

103.21 

51 

3. 7721 

111.35 

19 

C. 5996 

102 .33 

52 

0.95 50 

104.31 

20 

0. 1279 

99.79 

53 

0. 78 15 

103.63 

21 

1. 1196 

1C4 .90 

54 

0. 3489 

101.49 

22 

0.C368 

98.58 

3 5 

3.86 32 

111.52 

23 

0. 68Sd 

103.24 

36 

7.8031 ^ 

117.56 

24 

0.9047 

104.12 

57 

3. 26 3 4 

110.36 

25 

0.2734 

1C 1 .03 

58 

0.0253 

98.34 

26 

G. 87 15 

90.37 

59 

0.0146 

96.30 

27 

0.70 36 

iC3 .30 

60 

4.2446 

82.15 

28 

1.8903 

8 7,15 

6 l 

3.9288 

82 . 72 

29 

2.6336 

109.02 

62 

3. 2739 

8 3.98 

30 

0.00 16 

96.89 

63 

0. 71 31 

91.02 

31 

0. 5893 

91 .58 

64 

0. 1341 

94.05 

32 

2. 72 38 

35.14 

65 

0.94 73 

90.08 

33 

1 . 90 1 7 

37.12 

66 

0.^621 

102.14 


"S' OUTLIER IDENTIFIED AT 0.01 SIGNIFICANCE LEVEL 
SAMPLE SIZE IS 6 6 
NC OF OUTL lEP S IS 1 
DELSTAR = 6.3530 



DATA BEFORE IDENTIFICATION 

DATA AFTER DELETION 


OF OUTL lERS 

OF 

cltliers 

VARIABLES 

MEAN S .0 . 

MEAN 

S. D. 

H/ R 

97. 1BC6 7.2 956 

96. 8671 

6.8397 
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LIST OF VECTOPS WITH OUTLIERS IDENTIFIED BY ASTERISKS 


J 

DELTA 

G 

J 

DELTA 

G 

1 

1.3790 

0.70 

34 

0.8358 

0.85 

2 

1. 1829 

0.75 

35 

0.6848 

0.9C 

3 

1.0018 

0.80 

36 

0.5488 

0.95 

4 

0.8358 

0.85 

37 

0.4279 

1.00 

5 

0.68^8 

0.90 

38 

0-3220 

1.05 

6 

0.5488 

0.95 

39 

0.23 11 

1.10 

7 

0.4279 

1.00 

40 

0. 1553 

1.15 

8 

0.3220 

1.05 

41 

0.0945 

1.20 

9 

0.2311 

1.10 

42 

0.0487 

1.25 

10 

0. 1553 

1-15 

43 

0-0180 

1 .30 

11 

0.0945 

1.20 

44 

0.0022 

1 .35 

12 

0.0487 

1.25 

45 

0.00 16 

1 .40 

13 

0.0180 

1 .30 

46 

0 . 0 1 59 

1 .45 

14 

0.0022 

1.35 

47 

0.0897 

1 .55 

15 

0.00 16 

1.40 

48 

0. 1491 

1.60 

16 

0.0159 

1.45 

49 

0. 2236 

1 .65 

17 

0.0453 

1.50 

50 

0.3131 

1.70 

18 

0.0897 

1.55 

51 

0-4177 

1-75 

19 

C. 1491 

1.60 

52 

0-5372 

1.80 

20 

0.3131 

1.70 

53 

0.9861 

1.95 

21 

0.5372 

1.80 

54 

1.7951 

2.15 

22 

0.67 18 

1.85 

55 

2.03^9 

2.20 

23 

0.8215 

1 .90 

56 

4. 1346 

2.55 

24 

0.9861 

1 .95 

57 

4.4947 

2.60 

25 

1. 5703 

2 .10 

58 

1 4.1922 

3.55 

26 

1.7951 

2-15 

59 

1.00 18 

0.80 

27 

2.03^9 

2.20 

60 

0.8358 

0.85 

28 

2. 3137 

C .50 

6 1 

0.6848 

0.9C 

29 

2.0575 

0.55 

62 

0.5488 

0 .95 

30 

1. 5901 

0.65 

63 

0.4279 

1 .00 

31 

1.3790 

0.70 

6 4 

0- 3220 

1 .05 

32 

1. 1829 

0.75 

65 

0.0897 

1.55 

33 

1.00 18 

0.80 

66 

2-8446 

2.35 


* OUTLIER IDENTIFIED AT 0.01 SIGNIFICANCE LEVEL 
SAMPLE SIZE IS 66 
NC OF OUTL I ER S IS 1 
DELSTAR = 6.3530 



DATA BEFORE 
OF 

IDENTIF IC ATION 
OUTL lER S 

DATA 

AFTER DELETION 
OF OUTLIERS 

VARIABLES 

MEAN 

S .D . 

MEAN 

S 

G 

1.3773 

0 .576 7 

1.3438 0 
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LIST OF VECTORS WITh OUTLIERS IDENTIFIED BY ASTERISKS 


J 

DELTA 

H/R 

G 

J 

DELTA 

h/R 

G 

1 

1.4132 

91.12 

0.70 

34 

1. 2653 

89.42 

0.85 

2 

i.5222 

96.05 

0.75 

35 

0.9751 

9C.48 

0.90 

3 

1.3226 

96.32 

C.80 

36 

2.8389 

85.05 

0.95 

4 

3. 1165 

102.30 

0.85 

37 

0.5015 

96.03 

1.00 

5 

0.7199 

92.57 

0.90 

38 

0.6559 

98 .22 

1.05 

6 

0.6203 

95 .64 

0.95 

39 

0.2377 

95.63 

1.10 

7 

0.6784 

91.44 

1.00 

40 

3.4797 

106.36 

1.15 

8 

0.4002 

93.12 

1.05 

41 

0.0979 

96.23 

1.20 

9 

0.2635 

94.08 

I.IC 

42 

0.9158 

90.71 

1.25 

10 

0.5249 

91.90 

1.15 

43 

0.2202 

93.94 

1 . 30 

11 

0. 1091 

96.60 

1-20 

44 

0.3527 

93.46 

1.35 

12 

0. 1499 

98.14 

1.25 

45 

0-4353 

93.43 

1.40 

13 

0.6519 

101.35 

1.30 

46 

0. 1363 

95.65 

1.45 

14 

0. 1883 

99.55 

1.35 

47 

0. 19 67 

96.50 

1.55 

15 

0.0746 

95.74 

1 -4C 

48 

0.4257 

101.94 

1 .60 

16 

0.0263 

98.32 

1 .45 

49 

2.4135 

107.98 

1.65 

17 

0. 1421 

99.93 

1.50 

50 

0.3134 

99.45 

1 .7C 

18 

0.7313 

103.21 

1.55 

51 

4,1166 

111.35 

1.75 

19 

0.6055 

102.83 

1.60 

52 

0.9969 

104.31 

1.80 

20 

0.31 48 

99 .79 

1.70 

53 

1, 1297 

103.63 

1 .95 

21 

1. 1413 

104 .90 

1 . 80 

54 

1.6467 

101.49 

2.15 

22 

0.7922 

98,58 

1 .85 

55 

3.9882 

111.52 

2.20 

23 

0.9621 

103.24 

1 .90 

56 

8.0630 

117.56 

2.55 

24 

1. 1990 

104-12 

1.95 

57 

4.9986 

110.36 

2.6C 

25 

1-6294 

101.03 

2.10 

58 

20.3487 ={= 

98.34 

3.55 

26 

6. 1920 

90.37 

2.15 

59 

1. 3188 

96.30 

0.80 

27 

2.0352 

1C3.30 

2.20 

60 

4.36 16 

82.15 

0.85 

28 

2.6309 

67.15 

0.50 

61 

4.0832 

82.72 

0,90 

29 

11.1150 * 

109.02 

0.55 

62 

3.4158 

83.98 

0.95 

30 

2. 30 74 

96.89 

0.65 

63 

0. 75^1 

91 .02 

1.00 

31 

1.3906 

91.58 

0. 70 

6 4 

0.3372 

94.05 

1.05 

32 

2.7500 

•35.14 

0.75 

65 

2.0682 

90.08 

1.55 

33 

1.9632 

87.12 

0.30 

66 

2.9773 

102.14 

2.35 


OUTLIER IDENTIFIED AT 0.01 SIGNIFICANCE LEVEL 
S AMPL E S IZE IS 66 
NC OF OUTL IERS is 2 
DELST AR = 6.7115 



DATA BEFORE 
OF 

IDENTIFICATION 
OUTL IERS 

DATA AFTER DELETION 
OF CLTLIERS 

VARI ABLES 

MEAN 

S .D . 

MEAN 

S. D. 

H/ R 

97. 1806 

7.2556 

96. 9775 

7.2544 

G 

1.3773 

0.5767 

i. 3562 

0.5069 
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